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Abstract 

In this paper we will turn our attention to the problem of obtaining 
phase-space probability density functions. We will show that it is possible 
to obtain functions which assume only positive values over all its domain 
of definition. 
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1 Introduction 



Since the very moment, when the phase-space probability density representation 
was introduced to reproduce, in a manner as close as possible to classical sta- 
tistical mechanics, the results of the expectation values predicted by quantum 
mechanics |Q, there have appeared the discomfort of being capable of making 
this description only appealing to probability densities having negative values 
within some range of its phase-space arguments (x,p). 

Many papers have appeared since then trying to push this representation 
method to its utmost formal developments 0. The negativity problem was con- 
sidered by some as deriving from the non-commutativity character of operators 
corresponding to canonical conjugate functions, or just representing the fact 
that quantum mechanics, being an essentially non-classical theory, wouldn't be 
adjustable to a phase-space description. Others have thought this problem as 
irrelevant, for the phase-space distribution obtained (by Wigner-Moyal-Weyl 
procedure, for example) will furnish all the relevant expectation values identical 
with those of usual quantum mechanics. Indeed it will give even more, since 
expectations related with products of the type x n p m will be unambiguous cal- 
culated within this formulation of the problem — something the usual approach 
may do only after postulating some specific operator construction method. 

Others, however, have seen the far reaching consequences of being capable 
to construct such a representation, without the inconveniences introduced by 
the negativity problem || [|, 

We will join this last class of authors on their concerns. We have been 
arguing, since the first paper of this series^), that quantum mechanics is just 
a classical statistical theory, defined upon configuration space, related with the 
thermodynamic equilibrium situations of physical systems; this means that it 
is more restricted than the theory related with Liouville's equation, which is 
classical statistical mechanics strictu sensu. However, the element of liaison 
between the Liouville classical statistical mechanics, defined upon phase-space, 
and the Schrodinger quantum mechanics upon configuration-space was, until 
presently, the phase-space probability density function assuming negative values. 
If this function does assume negative values within some ranges of x and p values, 
then no acceptable physical interpretation may be appended to it. Then comes 
the discomfort of having all the theory based upon an element of irrationality 
(unphysical). 

The task of being capable to obtain positive definite probability densities will 
thus fulfill this epistemological void of the theory, giving to it its final contours. 
This is the intention of the present paper. 

In the second section we will show the method of deriving quantum me- 
chanics, we had thus far adopted, from classical statistical mechanics ||. While 
doing this, we will show what were our assumptions and why they did lead to 
the negativity problem. This critic development of the formalism will give us 
insight about the ways we may possibly overcome the difficulties. 
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In the third section we will use the insights furnished by the previous section 
to derive a method that gives us only positive definite phase-space probability 
density functions. 

The results of the third section will be used to derive, as a naive application, 
the Ehrenfest's theorem, giving the dynamical behavior of the position and 
momentum expectations. This will be done in the fourth section. 

In the fifth section we will compare the two procedures of finding phase- 
space probability density functions. This comparison will be made based upon 
the expectation values these two approaches furnish to the specific example of 
the harmonic oscillator. 

The last section will be devoted to our conclusions. 

In the appendix we make a comment about an interesting result, closely 
related with our own, made in the literature. 



2 Previous Approach 

In our previous approach to the problem, we began with the phase-space joint 
probability density function F(x,p;t) defined as a convolution of two phase- 
space probability amplitudes cj>(x,p;t) 

r + co 

F(x,p;t) = <f>\x,2 P - P ';t)<t>(x,p';t)dp'. (1) 



We then defined the characteristic function Zq(x, 5x;t) as the infinitesimal 
Fourier transform 

Z Q (x,Sx;t) = [ e^ Sx / h F(x,p;t)dp, (2) 



and applied the above defined Fourier transform to Liouville's equation, which 
we assumed the function F(x 1 p; t) satisfies, to find the equation satisfied by Zq 
as 

dZ Q (x,Sx;t) h 2 d 2 Z Q (x 1 8x\t) dV{x) 
~ lh at m dxd(Sx) +-^SxZ Q (x,Sx;t)=0. (3) 

Because of expressions (|l|) and |]), the characteristic function was given by 
Z Q (x,Sx;t) = ^ {x - 5x/2;t)*fj(x + 6x/2;t) (4) 
where we wrote, for the complex amplitudes, 

ip(x;t) =R(x;t)e lS{x '' t)/h , (5) 

with R(x; t) and S(x; t) real functions. The characteristic function became, in 
terms of the R(x; t) and S(x; t) functions 

Z Q (x,6x;t) = \R(x-t) 2 + 



3 



R{x;t) 



d 2 R{x;t) 
dx 2 



fdR{x;t) 
V dx 



exp 



i dS(x;ty 
h dx 



(6) 



up to second order in the variable dx, since it is assumed as infinitesimal and 
equation (||) has to be considered up to order one on this variable. 

Substituting these expressions in equation (|^) we got the two equations 



dR{x;tf 



dt 



d_ 

dx 



R(x;t) 2 dS{x;t) 



■m 



dx 



= 



(7) 



and 



-iSx d j dS(x;t) 



h dx 



dt 



J_ f dS{x;t) 
2m \ dx 



+ V(x) - 



d 2 R{x;t) 



2mR(x; t) dx 2 



= 0, 
(8) 

where we used the independence between x and Sx, and also the infinitesimal 
character of Sx to retain only up to power one on this variable. 

These two equations are mathematically equivalent to the Schrodinger equa- 
tion 



r d 2 ^(x;t) 



V(x)i>(x;t) = ih 



dij){x) t) 



(9) 



2m dx 2 dt 
when we substitute the expression (|^) on it, and separate the real and imaginary 
parts. 

This process allowed us to find, upon inversion of the Fourier transformation 
(^|), a phase-space probability density which will simulate all the results of ordi- 
nary quantum mechanics (the statistical moments, such as expectations, square 
deviations, etc.). However, by the very definition of the phase-space probability 
density F(x,p; t) given in ([!]), we shall not expect to find positive-definite func- 
tions, since the phase-space amplitudes may oscillate — mainly for the excited 
states, where their corresponding configuration-space amplitudes show such an 
oscillatory behavior — and their product is defined at different momentum-space 
points in (|l|), as required for a convolution. 

The process of operator formation was then simply introduced. Since Zq is 
a characteristic function (in momentum space) of the phase-space probability 
density function, all statistical moments may be obtained as 



(x n p ri 



= I lim x r ' 



—ih 



d 



8(8x) 



Z(x, 5x; t)dx, 



(10) 



where no ambiguity appears, since x and Sx are independent variables. Using 
the decomposition (0) we may connect these expectation values, calculated with 
the characteristic function Zq, with those calculated using the probability am- 
plitudes. This can be accomplished by means of the expression ([l(]) above with 
(||) for the characteristic function, after performing an expansion of (@) with re- 
spect to the infinitesimal variable Sx. In this case, the position and momentum 
operators become x = x and p = ihd/dx, and no longer commute. 
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We have argued that it was the process of inversion of the infinitesimal 
Fourier transformation (0) which introduced such negativity into our phase- 
space distributions. Now, we may say that it is mainly the imposition of having 
the characteristic function Zq written as a product that has introduced this 
negativity, since it forces the phase-space probability density to be written as 
a convolution in momentum space. We note, however, that, except for this 
unacceptable behavior — which has no physical interpretation — the phase-space 
densities thus derived was really adequate to the problem, for the reasons above 
mentioned. We just note that, when writing the characteristic function as the 
Fourier integral in (J^J) and assuming the decomposition in (||) we have also 
introduced an specific operator construction method which will give us a 
procedure to evaluating expectation values of functions as x n p m . This means 
that, when changing some of these definitions (e.g. the convolution in (|l|)), 
the operator construction method may also vary, giving possibly new expecta- 
tion values for these product. However, the expectation values of the 'lateral 
functions' x" and p m will always have to agree with those obtained by usual 
quantum mechanics techniques. 

We will now present a method to overcome the negativity problem. In the 
next section we will present a new derivation of the problem in which we assume, 
from the very beginning the positive character of the phase-space probability 
density function. 



3 Positive Densities 

We begin by defining the phase-space probability density as 

F{x,p;t) = $(x,p;t)(j>(x,p)t), (11) 

where <fi(x,p; t) are probability amplitudes defined upon phase-space. This func- 
tion will thus be obviously positive definite. 

The characteristic function continues to be defined by expression (|2|), but 
now we will also define the characteristic amplitudes 

£tf,x;t)= J e^'/ h 0(x,p;t)dp, (12) 

to write Zq, in terms of these functions, as 

Zq(x, Sx;t) = J £*(x', x; t)£(x' + Sx, x; t)dx', (13) 

which now is a convolution. This was readily expected, since now we are defining 
the probability density F(x,p;t) as a product upon phase-space, and F(x,p;t) 
and Zq(x, Sx;t) are the Fourier transform of each other. We thus have fixed 
the positivity of the phase-space probability density function while letting the 
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characteristic function to possibly assume negative values, since this will not be 
physically unacceptable. 

Another way of writing expression (13) is 

Z Q (x,5x;t) = J e{x',x;ty Sxp ' /h £,(x',x;t)dx\ (14) 

where 

p> = -m^. (is) 

We may substitute expression ( |l4| ) into equation (||) to find the equation 
satisfied by the characteristic amplitudes £ as 

- ffl^lM _ ^(x\x;t) dVix^ 

dt m dxdx' + dx x ^ x ' x ' l > u ' 

which was expected, since the phase-space amplitudes also satisfy a Liouville 
equation^ and the definition of the £'s parallels that of the characteristic func- 
tion Zq. 

It is obvious that all the statistical moments related to momentum or position 
values are still obtainable from the characteristic function Zq. In fact any 
function of the type 0[x,p] will have its expectation value, when calculated 
using the characteristic amplitudes, given by 



0[x,p] = I I C(x',x;t)0[x,p']£{x',x;t)dx'dx, (17) 

where 



d[x,p'] = o 



(18) 



and no ambiguities appear, for x and x' are independent variables. These results 
come from the very application of the definitions ( |l0| ) to the new definition ( |l3| ) . 
The explicit results are 

{x n p m ) = J j \{x' ,x;t)x n f-ih^\ Z(x',x;t)dxdx'. (19) 

We may now put 

£(x', x;t) = i){x' \t)ijS {x;t) (20) 
to get, for the characteristic function 

Zq(x, 5x; t) = ip(x; t)ip^ (x; t) J ij^(x'; t)exp ^Sx— tp(x'; t)dx' , (21) 

giving, for the lateral statistical moments 

W) = J ^{x'\t){-i%—\ i>(x';t)dx' (22) 



G 



and 

J^)= I ^{x)t){x n )ij]{x;t)dx, (23) 



whereas, for the general product, we have 

{x n P m ) = J ^(x]t)(x n )ip(x;t)dx J ^(x';t) (~ ih Tj^) i>(x';t)dx'. (24) 

Note that, now, the characteristic function became a quadrilinear function of 
the amplitudes^. This new definition of the characteristic function will lead 
to an operator construction method distinct from the one presented in the last 
section. 

Until now we have not specified which equation the probability amplitudes 
ip(x]t) satisfy — indeed, this was done in the previous section but, since our 
definitions has changed, we are not certain if these amplitudes still obey the 
Schrodinger equation. However, if we turn our attention to the operator forma- 



tion procedure given in (|22j) and (23), we thus see that the expectation of the 
energy function, given by 

E = J J (j- + V(x)jF(x,p;t)dxdp, (25) 

may be written, in terms of the probability amplitudes, as 



E = J i>\x;t) 
which implies the Schrodinger equation 



1 / d ^ 



rp(x;t)dx, (26) 



j^iP(x,t) + V{x)iP{x;t)=Ei>(x;t), (27) 

with p — —ihd/dx. This suffices to prove that the amplitudes ij)(x;t) still obey 
the Schrodinger equation. 

The last expression for the expectation values of coordinate-momentum pow- 
ers ( p^ ) will thus strongly differ from the one we have obtained in the last sec- 
tion. We also note that, in the present procedure, we still do not introduce 
any commutation problem, for x and x' continue to be independent variables, 
just like x and Sx are. Since we showed that the probability amplitudes sat- 
isfy the Schrodinger equation, all the 'lateral' expectation values predicted by 
quantum mechanics will be reproduced by the present theory. The discrepancy 
between this section approach and the one described in the previous section 
may not be solved within the realm of the usual quantum formalism, for both 
approaches predict the same lateral statistical expectations (x n or p m alone), 
which are the only ones guaranteed by the above mentioned formalism — the 
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coordinate-momentum product expectations depending on the particular oper- 
ator construction method (the correspondence rule) each approach proposes. 
Using now the definition (\L4\) an( l the decomposition (EQ), we have 



^(x;t)ijj(x';t) = e lpx ' /h <j)(x,p;t)dp, (28) 



<l>(x,p;t)=^(x;t)ip(p;t), (29) 



(f(p;t)= / e- ipx '' n i>{x';t)dx'. (30) 



which implies 
where 



With the result (29|), for the phase-space probability amplitudes, and the 
initial imposition (11) we get 

F n {x,p;t) = \i> n {x;t)\ 2 \y n {p;t)\ 2 , (31) 

as the desired positive definite phase-space probability density function. This 
result reinforces the F n property of giving all the lateral statistical expectations 
as the ones predicted by usual quantum mechanics. 

In the literature there are some restrictions to the density obtained in (|3l|). 
The nature of these restrictions^ and the explanation of why we have been 
successful in overcome them with our present approach may be found in the 
appendix. 



4 Ehrenfest's Theorem 



From the above relations it is possible to derive the Ehrenfest's theorem in a 
straightforward manner. To see this consider first the equation satisfied by the 
characteristic function Zq (^) . Writing this characteristic function up to second 
order in the infinitesimal displacement 5x we get 



Zq(x, 8x\ t) — ijr{x\ t)ip(x; t) 



(32) 



Now, taking this expression into (|3j) and separating the zeroth and first order 
terms in the infinitesimal displacement 8x, we get the two equations 



in^(x;t)^(x;t)) - ^-^-(^(x;t)^(x;t)) f ^ (x'; ^^ 



dx' = 0, 
(33) 
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and 



-ih{i)\x\t)il){x;t)) J ip{x';t) 
h ~ d -{^{x;t)i>(x;t)) f ^(x' ^t^ dx' 



dx> 
di/j(x';t) 



dx 



dx' 2 



dx 

dV(x) 
dx 



Multiplying the first equation by x and integrating we find 

dx 



Of 



p. 



while, integrating the second equation in x, we find 



dp 
di 



dV(x) 
dx 



dx' — 

{^{x;t)i>(x;t))=Q. 

(34) 

(35) 
(36) 



These last two equations are exactly the mathematical expressions of the Ehren- 
fest's theorem]?]]. 

Moreover, with the probability density on configuration space written as 
p{x;t) — \ip(x;t)\ 2 , the equation ( |33|) becomes 



dp(x;t) d_ 
dt dx 



P(t) 



p{x;t) 



0, 



while using this equation into (|34| ) we get 
'dpltj dV(x) 



p(x;t) 



dt 



dx 



(8pfdp(x-i) 
m ox 



(37) 
(38) 

111 VJu 

JW = pW~pJt) 2 - (39) 

These last two equations are precisely those we would get if we integrate the 
Liouville equation in the variable p after multiplying it by 1 and p, respectively, 
using a separable density function of the type given in (|l]) (l4| . 



where 



5 Application: Harmonic Oscillator 

We now want to make a simple application of the previous results to the har- 
monic oscillator problem. Thus, we begin with the harmonic oscillator proba- 
bility amplitude given by 

M*',t) = (a 2 /4) V2 ^-^'\-^l 2 U n{0Lx)e ^l\ (40) 



where H n (x) are the Hermite polynomials and a — y mco/h, with m the mass 
and to the frequency related with the problem. 

Using the approach of the second section, the phase-space probability density 
function may be written as[^| [l0| O, [l2| 

F*(x,p;t) = ^(-l) n eay (~\ L n (^j , (41) 

which clearly assumes, for the excited states, negative values within some ranges 
of the variable H = p 2 /2m + muj 2 x 2 /2. 

The phase-space probability density function of the present approach we 
uphold in the third section is simply 

F*{x,p;t) = \M^t)\ 2 \Vn{p;t)\ 2 , (42) 

where cp n (p; t) is the solution of the problem on the momentum space — the 
Fourier transform of i/j(x; t). 

The results obtained for the expectation values by these two phase-space 
probability densities are exhibited in Table I. We may see from this table that 
the expectations will coincide for all cases but those involving the product x n p m . 
While the phase-space probability density function (^) presents some correla- 
tion between the variables x and p, these correlations are absent from the one 
in ([l2|). In the present harmonic oscillator case, this will affect, for example, the 
expectations of the energy dispersions, which become higher for the uncorrelated 
density. 

This result was expected since, for the excited states, the densities spread 
more and more with respect to the variable H and its dispersion, calculated with 
the probability density F B , remains constant only because this density assumes 
negative values within some ranges of H; since this negativity regions do not 
appear in the probability density F A , we may expect the spread in the energy to 
get larger, following the density spread itself. Thus, the probability density F B 
presents the strange behavior of having the momentum and position dispersion 
getting larger and larger, while the energy dispersion being kept constant. 

We have also plotted the first three phase-space probability density functions 
according to the prescription (p2[). They are given by figures I, II and III. These 
figures may be contrasted with those related with the prescription ( fi"l| ) (see, for 
example, QlavoJl3|). 

6 conclusion 

In this paper we have been successful in obtaining positive definite phase-space 
probability density functions. Although this result does not change very deeply 
the formal aspects of the quantum theory, since this theory will continue to be 
given by the Schrodinger equation, it does constitute an important epistemo- 
logical acquisition. Indeed, since the very beginning of this series we have being 
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arguing that quantum mechanics is a mere classical statistical mechanics, per- 
formed upon configuration space, for systems in thermodynamic equilibrium E3 . 
However, it will be rather inconvenient for such an epistemology to have its main 
symbol of liaison (the phase-space probability density) between these two lev- 
els of descriptions (phase-space and configuration space) not being capable of 
having a physical interpretation. 

Thus, being capable of calculating the probability density functions on phase- 
space not suffering from the negativity problem gives us the final justification to 
our previous assertions. Now the theory has all its symbols with some acceptable 
physical interpretation attached to it and we may finally say, without fearness, 
that quantum mechanics is just an ensemble statistical theory performed upon 
configuration space and related with thermodynamic equilibrium situations. 

The next paper will be devoted to the epistemological consequences of the 
formalism derived in this series of papers || and will constitute its natural logical 
closure. 



A Restriction upon the Density 

Cohen || has derived a generalized method of obtaining phase-space density 
functions. He was aiming at giving a general procedure of obtaining all the 
phase-space probability distributions proposed at the time his paper was pub- 
lished. He thus begins by writing this phase-space probability density function 
as 

F(x,p;t;f) = ^ J J J e- i6x - iT P +ieu f(6,r;t)x 

x^(u~ ^rh; t)ip(u + -rh) t)d6drdu. (43) 

In this case, a wide variety of probability densities may be obtained by just 
fixing the function f(9,r;t) (e.g. the Wigner function in (|^) may be found by 
fixing / = 1). 

In his conclusions he argues that: "It is commonly held that the uncertainty 
principle by itself precludes the possibility of the existence of a joint distribution 
of position and momentum. However, this is not so. For example, the choice 

= f f\^(x)m P )\ 2 e^ Pdxdp 

J( ' T ' ' f^(u-±Th)e^(u+±Th)du ( ' 

leads to 

F(x,p;t) = |^(x)| 2 |0(p)| 2 , (45) 

which is certainly a well-defined joint distribution and from which the uncer- 
tainty principle follows in the usual manner. The reason why a true joint distri- 
bution cannot be defined is because no choice of / yields a distribution which 
gives the correct quantum mechanical expectation values for all observables 



11 



when calculated through phase-space integration. That is, no / exists such 
that, if the correspondence of quantum to classical variables (..) is 

g(x,p)^G, (46) 

for some /, then also 

H(g(x,p)) - H(G), (47) 

for the same /, where H is any function" . 

However, this was exactly what we have obtained in the third section of this 
paper. Indeed, our method of operator formation was 

d 

x = x and p = —ih— — (48) 
ox' 

for which any function H will be such that 

H(g(x,p)) - H{G) = H(g(x,p)), (49) 

since x and p commute, denying ambiguities related with operator ordering to 
appear. 

This is precisely the point of disagreement. When constructing his correspon- 
dence rule, Cohen[[l5) has explicitly assumed that the operators related with the 
position and the momentum are those given by x = x and p — —ihd/dx, which 
do not commute. This precludes, in his approach, the possibility of finding a 
function / which gives, at the same time, a specific correspondence rule, the 
density as in ( [45| ) (for example), and all the expectations values predicted by 
quantum mechanics as equal to those calculated by using the derived correspon- 
dence rule. 

Thus, with the choice ([l3|) he was able to find a function / which gives 
the phase-space probability density as ([45|). With this function / he was able 
to find the correspondence rule which takes a function g{x,p) into its related 
operator G. But when calculating the expectation values with the phase-space 
probability density and the operator G, he was not able to make the results 
agree. Looking at our own development, we may see that the problem is due to 
the fact that, in his last derivation step, Cohenjl5| imposes that the variables x 
and p be substituted by the non-commuting operators x — x and p — —ihd/dx, 
which may not be accommodated within our present approach. Indeed, it is 
the commuting of the operators in (^8|) that guarantees the validity of the ex- 
pectation values calculated with the uncorrelated probability density (^5|), as 
we can see by (p4|), and it is obvious that we will have exactly the same results 
while calculating the expectations with the phase-space probability density F A 
or using the operator given by (f49|), which is our correspondence rule. 
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1/4 


3 
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7/2 


7/2 


75/4 


49/4 


25/4 


13/4 


1/4 



Table 1: Some expectation values calculated with the present probability density 
function in phase space (A = F A ) and the one usually found in the literature 
(B = F B ). The results coincide whenever the index is absent. 

TABLES 
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FIGURES 



Figure 1 - Plot of the phase space probability density function for N=0. 



Figure 2 - Plot of the phase space probability density function for N=l. 



Figure 3 - Plot of the phase space probability density function for N=2. 
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FIGURE I 




FIGURE II 



FIGURE III 




